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Abstract. Cygnus X-3 is a strong X-ray source (Lx ~ 10^^ erg s~^) which is thought to consist of a compact object, accreting 
matter from a helium star. We find analytically that the estimated ranges of mass-loss rate and orbital-period derivative for 
Cyg X-3 are consistent with two models: i) the system is detached and the mass loss from the system comes from the stellar wind 
of a massive helium star, of which only a fraction that allows for the observed X-ray luminosity is accreted, or ii) the system 
is semidetached and a Roche-lobe-overflowing low- or moderate-mass helium donor transfers mass to the compact object, 
followed by ejection of its excess over the Eddington rate from the system. These analytical results appear to be consistent with 
evolutionary calculations. By means of population synthesis we find that currently in the Galaxy there may exist ~ 1 X-ray 
binary with a black hole that accretes from a > 7Mq Wolf-Rayet star and ~ 1 X-ray binary in which a neutron star accretes 
matter from a Roche-lobe-overflowing helium star with mass < 1.5Mq. Cyg X-3 is probably one of these systems. 
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1. Introduction 

Cygnus X-3 was discovered as an X-ray source by 
tGiac coni et all It is a strong X-ray source (Lx ~ 10 

i erg s~^, assuming a distance of 9 kpc), and the X-rays are 
' expected to be due to accretion of matter onto a compact 
object (c.o.), presumably a black hole ( BH) or a neu tr on sta r 
(NS) (see, e.g., iKitamoto et'aLl Il987t IPredehl et all \2mi . 
The X-ray and infrared (IR) emission show a periodicity of 
4.8 hours, which is believed to be the orbital period P of the 
system (see, e.g., IParsignault et alJ Il972h . Van den Heuvel 
& de Loore (1973) suggested that Cyg X-3 consists of a NS 
with a helium (He) star companion, as a later evo lutionary 

froduct of a high-mass X-ray binary. iTutukov & YungelsonI 
1973ah independently considered a NS accompanied by a 
He star as a stage in an evolutionary scenario leading from 
a pair of main-sequence stars to a binary NS. There is too 
much interstellar obscuration towards the source to observe 
it optically, but observations in the IR wave bands in the 
1990's by van Kerkwijk and coauthors (1992, 1993, 1996) 
identified a Wolf-Rayet (WR) spectrum with Cyg X-3. Both 
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the observations of van Kerkwijk and coauthors as well as 
high-resolution spectroscopy by Fender etal. (1999) revealed 
hydrogen depletion of the mass donor. Furthermore, phase- 
to-phase variations in the X-ray spectra can be explained by 
a strong (fact or 10-100) overabundance of carbon, nitrogen, 
or oxygen ( iTerasawa & Naka mura 1994), consistent with a 
classification of the Cyg-X-3 companion a s a WN-typ e star. 
This added credibility to |van den He uvel & de Loora s and 
iTutukov & YungelsorJ s prediction. 

The aims of this paper are twofold. 

(1) We will determine what combinations of stellar masses 
and what mass-transfer/mass-loss mechanisms are consistent 
with the observed P /P values and observed mass-loss rates 
M from the Cyg-X-3 system (the observations are taken from 
the literature and summarised in Sect. |2j- For this we make 
analytical estimates and carry out evolutionary calculations for 
systems consisting of a He star and a c.o. (Sects. |3and|3}. 

(2) We will study how many X-ray sources with these (Cyg- 
X-3-like) parameters are expected to exist in the Galaxy at 
present. To do this, we will caiTy out a population synthesis for 
He-star plus c.o. (He-nc.o.) binaries in the Galaxy and apply 
disk-formation criteria to estimate the number of observable 
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X-ray binaries with Cyg-X-3-Hke parameters in the Galaxy 
(Sect. |5|l- We try to explain why Cyg X-3 is the only X-ray 
binary with a He-star companion we observe in our galaxy 
(Sect. 16. 1> , and briefly discuss two other binary systems that 
were recently suggested to consist of a c.o. and a WR star 
(Sect. l6T2t . A Conclusion follows in Sect.0 

2. Observed and inferred properties of 
Cyg n us X-3 

2. 1. Masses of the components 

iTavani et alJ ( Il989l) and later iMitral ( Il99d Il998l) argued that 
the companion in Cyg X-3 should be a very-low-mass (0.01 — 
O.OSMq) He star. However, this is hard to reconcile with the 
high IR luminosity of the system ( Hanson et al. 200 0). Our cal- 
culations (Sect. l3.2l &l4ll and population synthesis (Sect.|5} also 
exclude a donor less massive than 0.8Mq. 

[Ter asawa & Nakamura (1994) found, from the ionisation 
structure of the wind in Cyg X-3, that the mass of the wind- 
supplying component has to be moderate: Tlj^^©. 

ISchmutz et al. (19963 conclude that the variations in the 
profiles of several near-IR emission lines are due to the or- 
bital motion of the WR star and derive a mass function for 
the donor f{m^) — 2.3 Mq. For the range of assumed Wolf- 
Rayet masses 5 to 20 Mq and a range of possible inclinations 
30° < i < 90° they get a mass in the range (7 — AO^AIq for 
the CO. from which they infer that it is a BH. 

iHanson et alJ ( l200(]l) found a sinusoidal absorption feature 
originating in the wind in the 2.06/im spectral region of Cyg X- 
3 which allowed them to derive a mass function /(m) = 
0.027 Mq. Hanson et al. considered two options: an origin of 
the absorption in the accretion disk or other material centered 
on the compact object, or association of absorption with the 
donor. The first option is consistent with low- or moderate- 
mass (< 8 Mq) donors, but requires a low orbital inclination of 
the system (< 20°). Association of the absorption feature with 
the donor limits the mass of the WR companion to < IOMq 
if the accretor is a NS. For BH accretors the mass of the sec- 
on dary may be a s high as 7OM0. 

IStark & Saial ll2003t) studied the modulation of X-ray emis- 
sion lines from Cyg X-3. Based on a discussion of the location 
of the regions of emission of highly-ionised silicon, sulfur, and 
iron, they assume that the iron line is produced in the wind 
captured by the c.o. or in an accretion disk around it. They then 
use the fact of non-detection of a modulation of the iron lines 
to derive an upper limit to the mass function for the accretor: 
/(ma) < O.22M0. For an accretor of 1.4 M© or 10 Mg, the 
minimum mass of the donor is then ~ 1.1 Mq and ~ 3.4 Mq, 
respectively. We furthermore note that in the case of Roche- 
lobe overflow (RLOF) a mass ratio nid/nii^ > 1.39 is incon- 
sistent with the observed increase of the period (see below), 
which then also implies an inclination > 20° in this case. 

The early models of Cyg X-3 that assumed an elliptic or- 
bit for the system (Ghosh et al. 198 1 ) may be discarded now, 
since no sign s of an apsidal mo tion were found in ~ 30 yr of 
observations ('Singh et alJl2002h . An implication of this non- 
discovery of an apsidal motion is the irrelevance of the values 



Table 1. The values for P/P of Cyg X-3, derived by fitting two 
different models to the observations. A parabolic ephemeris as- 
sumes P = 0, whereas in the cubic ephemeris also a second 
derivative of the period unequal to zero is taken into account. 



P/P, fitted value 


X?od Reference 


(10-® yr-i) 




Parabolic ephemeris 


2.2 ±0.3 ^ 


1.41 van der Klis & Bonnet-Bidauil ^1981) 


2.19 ±0.05 t 


0.78 Kitamotoetal. fl987) 


1.6 ±0.1 


1.55 van der Khs & Bonnet-Bidaud (1989) 


1.2 ±0.4 


2.07 Kitamoto et al. (1995) 


1.05 ±0.04 


3.08 Sineh et al. (2002) 


Cubic ephemeris 


4.0 ±0.6 1" 


1.36 van der Klis & Bonnet-Bidaud ( 1989) 


2.9 ±0.2 + 


1.46 Kitamotoetal. (1995) 


1.4 ±0.3 + 


2.96 Sinsh et al. (2002) 



f - value not quoted in the referred-to paper, but calculated by us from 
the published values of P and P. 



of the orbital inclination found bv iGhosh et alJ ( Il98lh . which 
are often used in the literature. 

To summarise, a great ambiguity still exists in the interpre- 
tation of the radial-velocity curves of Cyg X-3, mostly related 
to different locations of spectral features that serve as the basis 
for radial- velocit y determinations. H owever, at the moment it 
seems Ukely that lStark & Saial ( l2003l) really measure emission 
originating in the vicinity of the c.o. Then their results suggest 
a rather moderate mass for the companion to the c.o., if the c.o. 
is a neutron star or a stellar-mass black hole. 

2.2. Orbital period and its derivative 

The period P of Cyg X-3 has been e xtensively moni- 
tored over the years (' e.g., lvan der Klis & Bonn et-Bidaud 19*811 
iKitamoto et alJll987l) . and is found to be increasing on a rela- 
tively short time scale of ~ 10^ yr. There are also indications of 
a second derivative of the order of — 10~^° vr~^ to the per iod 
dvan der KUs & Bonnet-Bidaudll989HKitamoto et alJl995l) . A 
summary of the estimates of P/P for Cyg X-3 is presented in 
Tabled 

2.3. Mass-loss rate 

The mass-loss rate for Cyg X-3 was estimated from IR ob- 
servations, usually using the Wright & Barlow ("1975", W & B) 
model for the emission of a spherical, homogeneous, constant- 
velocity, isothermal wind. Stars have an accelerating wind with 
a temperature gradient, but W & B note that observations show 
spectrum flattenings in the near-IR similar to those predicted 
by their constant-tempera ture , constant- velocity model . 

Waltman et al. Vl996) and'Mill er-.Tones et alJll2005l) use an- 
other method to estimate M. They also assume the mass out- 
flow to be spherically symmetric and then use the fact that a 
(post-outburst) jet becomes observable with different delays af- 
ter the burst at different frequencies (see .Waltman et aL.19%1 
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Table 2. Values for M from the Hterature. We only show values 
that are obtained from observations of the mass loss from the 
system, hence not those inferred from the evo lution of the or- 
bital period. All estimates except lOglev et all (12001.) '^ assume 
spherical symmetry. See the main text for details. 



Estimated M 

(0.2 - 2.7) X IO-'^Mq yr"^ 
4 X lO-'^Mo yr-i 

< 10-*Mq yr-i 

(0.4 - 2.9) X IQ-^Mq yr-i 

< 10~^Mq yr^i 

~ 1.2 X 10"'' Mq yr"^ 
(0.5-3.6) X lO-^Moyr-i 



Reference 



Waltman et al. (1996) =^ 
van Kerkwiik (•1993) '^ 
van Kerkwiik et ak (1996) '' 
Oglev et al. (20^ 
Oelev etal. (2001) *^ 
Koch-Miramond et al. (2002)'= 
Miller-Jones et al. (2005)' 



^ From delays between the 15, 8.3, and 2.25 GHz radio light curves, 
assuming a jet velocity of 0.3 c. 

^ W & B model, from K-band observations, taking the wind veloc- 
ity to be w„ = 1000 km s"^ 

W & B model, from new I- and K-band observations which gave 
an improved value for the wind velocity ~ 1500 km s~^. 
Adojgting the non-s pherical, disk-like model by iFender et alj 
il992) and using the lGorensteinI il975h approximation for X-ray 
absorption. 

W & B model, taking = 1500 km s"^ 
' From delays between the 43 and 15 GHz radio light curves, as- 
suming a jet velocity of 0.6 c. 



for details). Note that this method gives M in the wind, not M 
in the jet itself, which is assumed to be much smaller 

The estimated mass-loss rates for Cyg X-3, varying from 
0.5 X 10~'^AfQ yr~^ up to 2.9 x IO^'^Mq yr"\ are presented in 
Table|2] All estimates except one assume spherical symmetry. 
Note that deviations from spherical symmetry will most proba- 
bly resi ilt in a higher mass-loss rate in the estimates f rom time 
delays ("Waltman et al. 1996: Mil ler-.Tones et al.l2005l) . whereas 
deviations from spherical symmetry in the other cases will re- 
sult in a lower effective mass-loss rate from the system (see, 

I » I 

e.g., Ko ch-Miramond e t al. 2002). 

If the increase of the orbital period is considered to be the 
result of a high- velocity wind from the system that takes away 
specific angular momentum of the donor (e.g..'Kitamot o et all 
1 1995: Ergma & Yungelson ,1998) the formula P/{2P) = 
-M/Mt yields M 5 x lO-'^(Mt/M0)MQ yr-\ where Mt 
is the total mass in the system. 



3. The models 

We consider two possible mechanisms that may cause varia- 
tions of the orbital period in a binary system consisting of a 
CO. and a companion. In the first model the companion loses 
mass in a wind which is directly lost from the system, ex- 
cept a tiny fraction that may be accreted (we assume accre- 
tion at the Edding ton limit). This is the mod el assumed by, 
e.g., iKitamoto e~al. ( 1995 ) and Erg ma & Yungelson ( 199 8), 
see previous Section. In the second model the companion trans- 
fers mass to the c.o., which then may eject (part of) the trans- 
ferred mass from the system. A similar quahtative model was 



first suggested for Cyg X-3 by IVanbeve ren et a l.l (Il998al) and 
later for both SS 433 and Cyg X-3 by Fuchs et alJ ( l2002al 
2004). In this model the transferred matter forms an envelope 
around the c.o., resembling a stellar atmosphere. A small thin 
accretion disk may be present around the c.o. The envelope is 
ionised by X-ray emission from the vicinity of the c.o. and ex- 
pelled from the system by radiation pressure. The whole pro- 
cess mimics the formation of a WR-star wind and consequently 
of a WR-like spectrum. Note that a disk-like wind from the sys- 
tem would give the same observed spectrum as a direct spheri- 
cal wind from the donor Following Fuchs et al. (2002a[ Eo04l) 
one may speak about a "WR phenomenon" in this case. The 
line-emission region may then really be associated with the 
c.o. as suggested by, e.g.. Fender et al. ( 1999) and Star k & Saial 
(2003). The details of this qualitative picture remain to be elab- 
orated upon and verified by observations. 

In principle the observed mass loss may be a combination 
of a direct wind and re-ejection of transfeiTed mass. However, 
since in the case of Roche-lobe overflow (RLOF) the mass-loss 
rate considerably exceeds the mass-loss rate of a direct wind 
(see Sect. 13. 2> . we consider only the extreme cases, one being 
only wind mass loss (with only so much mass transfer as to 
allow the c.o. to accrete at the Eddington rate), the other being 
no direct wind and only mass transfer from the He star followed 
by re-ejection from the compact star. 

3. 1. Equation for the period derivative 

We consider a binary system containing a mass-losing star (the 
donor) with mass md and a c.o. (the accretor) with mass rria- 
The donor loses an amount of mass drrid of which a fraction a 
is directly lost from the system, carrying away the specific or- 
bital angular momentum of the donor A fraction /3 is first trans- 
feiTed to the accretor and then lost through re-ejection, caiTying 
away /3dmd in mass with the specific orbital angular momen- 
tum of the accretor. We can then find P/P as a direct function 
of a, f3, rrid, nin, and riid (see Appendix A): 



md 



P 



X [(3m^ — 2maTOd — irin?^a — 2ma?nd/3 + 3md — 3m 



(1) 



We assume the c.o. to accrete at the Eddington rate, which 
leaves four free parameters to deal with to obtain a time scale 
for orbital-period variations in the range of observed values for 
P/P (see Table 0. Note that the cubic ephemeris (P ^ 0, 
in all those cases P < 0) fit the observations (slightly) better 
th an the parabolic ephemeris. Although the c ubic ephemeris 
of Ivan der Klis & Bonnet-Bidaud (1 1989!) and IKitamoto et"al] 
("1995^ are not confirmed — according to both ephemeris the 
period P should have stopped increasing and started decreas- 
ing by now, which is not observed — we note that P < may 
indicate a decrease in the mass-loss rate. Discarding the not- 
confirmed cubic models c and d in Table[2 the observed P/P 
values are in the range 1.0 x 10^® yr^^ to 2.2 x 10^® yr^^. 
We note here that in the extreme case in which there is only 
re-ejection and no direct WR wind (i.e., a — (3 — 1), only 
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Fig. 1. Mass-loss rates of He stars: (a) - wind mass-loss rate 
for a homogeneous He star; (b) - wind mass-loss rate for a He 
star at He terminal-age main sequence; (c) - rough estimate of 
the mass-loss rate of a star that overflows its Roche lobe on a 
thermal time scale after completion of core He burning. 

ratios of md/nia < 1.39 give an increase of the orbital period 
and are thus consistent with the observed P of Cyg X-3. 

3.2. Limits on tine component masses derived from 
ttie rate of mass transfer and mass loss 

3.2.1 . Mass loss due to a direct wind 

The presence of relativistic jets (e.g., iMioduszewski et all 
1200 1*) shows that there must be an accretion disk around Cyg 
X-3. This accretion disk may be either formed through wind ac- 
cretion or through RLOF. We will first consider the case where 
the He star loses mass in a wind and the compact object ac- 
cretes part of this wind (presumably at the Eddington rate). For 
the wind mass-loss rate of a He star we use 



_ r 2.8 X 10-i3(L/Lo)i-5 



if logL/Lo >4.5 

if log L/L© < 4.5 ^ ' 



(in Mq yr ^, see IPewi et aDl2002l and references therein). 
We take the lu minosity of a He sta r on the He m a in-seq uence 
(HeMS) from iHurlev et al.1 tOOd . iHurlev et al.1 (l200oi) only 
did c alculations for He stars up to l OMf^. but'Nugis & LamersI 
1I2OOQ1) find that when lHurlev et al.f s results are extrapolated up 
to He-star m asses of Mhp = 40Mf7> the y are consistent with 
the results of ISchaerer & MaedeJ (Il992'). which are valid for 
WNE'AVC stars up to 40Mq. In Fig. [T]we plot the thus ob- 
tained mass-loss rate for a He star as a function of its mass 
Mb_c at the beginning and at the end of the HeMS. 

Observed populati on I WR stars ar e more massive than 
~ 7Mq (see, e.g., Nug is & Lamersl2000l) : lower-mass He stars 
will probably not show a WR spectrum, nor will they pro- 
duce a wind that cou ld explain the ma ss-loss rate observed in 
the Cyg-X-3 system. iNugis & LamersI do not report any WR 
stars with mass-lo ss rates below 4 x 10~^A/q yr^^, whereas 
iMiller- Jones et alJ (.2005) find that the mass-loss rate from Cyg 

' early WN stars, WN2-5 



X-3 may be below 4 x lO^^M© yr^^. This might be an indi- 
cation that the WR spectrum observed from Cyg X-3 is due to 
the re-ejection of Roche-lobe-overflowed material that mimics 
the WR phenomenon. Another explanation, however, may be 
that the outflow is not spherically symmetric in the case of Cyg 
X-3 (Sect.|3. 



3.2.2. IVIass transfer due to RLOF 

Evolutionary calculations of 'Paczvriski' iIh 
ljutukov & Yungelson ( 19731^, and |lben & Tutuko\] lfl985l) 
showed that He stars with A/hc S O.SA/q do not expand 
during core He burning and later evolutionary stages. Also 
He stars with A/hc ^ 8Mq hardly expand before carbon (C) 
ignition in their cores (later stages are so short that they may be 
neglected). For intermediate-mass He stars, inspection of the 
summary figure of Dewi et al. (2002, their Fig. 1) shows that 
in a binary with P — 0.2 days containing a NS and a He star, 
the latter may overflow its Roche lobe in the He-shell-burning 
stage if A^Hc ^ 5.8Mq or in the core-C-burning stage if 
-^^Hc ^ 7.4Af0. However, the expected number of systems 
in the Galaxy that experience RLOF in the C-burning stage is 
negligibly small since this stage is very short, and we are left 
with (0.8 - 5.8)Afo He stars, overflowing their Roche lobe in 
the He-shell-burning stage (so-called BB case of evolution). 
If RLOF occurs after core He burning is completed, the 
mass-exchange time scale is of the order of the thermal one: 



GAP 



Tth = 



RL ' 



(3) 



with M, R, and L the mass, radius and luminosity of the donor. 
In Fig. [1] we also plot the mass-transfer rate from the He star 
given by Eq. (O, with R and L at the end of the He main- 
sequence from Hurlev et al. ( 20QQJ . We see that if in a system 
containing a rather low-mass He star matter is transferred to the 
CO. via RLOF in a thermal time scale, the rate of loss of the re- 
ejected matter from the system would cover the same range in 
Al as that for the winds of WR stars. Notice that Fig.fQjustifies 
the consideration of only the two extreme cases of mass loss 
(stellar wind vs. RLOF-nre-ejection): in a system with a given 
He-star mass that experiences RLOF, the mass transfer rate due 
to RLOF is two orders of magnitude larger than the wind mass- 
loss rate of a He star with the same mass would be. 

3.2.3. Analytical results for the Cyg-X-3 system 

We can now insert the estimated values of RLOF and wind 
mass-loss rates given by Eqs. (|2} and (|3} and the observed P/P 
into Eq. Q and solve it for all possible combinations of masses 
of the components rria and md that satisfy the observations for 
Cyg X-3. The result of this is shown in Fig.|2 The narrow strips 
in the lower left show the possible combinations, assuming that 
mass transfer in the system is due to RLOF in a thermal time 
scale. The strips are shown here to exemplify their narrowness: 
the width of the strip is typically less than COSAf©. When one 
uses mass-transfer rates for RLOF obtained from real binary- 
stellar-evolution calculations, the position of the strip may be 
somewhat displaced, but it remains narrow. The larger areas 



Dave Lommen et al.: Cyg X-3: the only WR X-ray binary 



5 



o 10 

03 







wind accretion 


llllllllllllllllllllllll' " 






































iNb systems 



5 10 15 

Donor mass (M©) 
Fig. 2. Schematic representation of possible system configu- 
rations that are consistent with P/P and M as observed in 
Cyg X-3. RLOF solutions give a very narrow range of pos- 
sible combinations of and toj (the two narrow strips in 
the lower left; the gap is due to the presumed gap between 
the masses of NSs and those of BHs). For wind-accretion so- 
lutions the range of possible combinations of and is 
much larger (hatched areas). Population synthesis shows that 
RLOF systems with NSs are formed in sufficiently large num- 
bers to produce an observable Cyg-X-3 system. On the other 
hand, Roche-lobe overflowing BH systems are formed so rarely 
with P similar to the period of Cyg X-3 that the probability of 
observing them is negligible. See Figs.lsjandl^and text. 

to the right show the possible mass combinations, assuming 
that mass-transfer in the system is due to wind accretion. We 
note that the much larger areas for the wind-accretion solutions 
do not mean that wind accretion is the most probable means 
of mass transfer in the Cyg-X-3 system. It merely shows that, 
;/ wind accretion is the mass-transfer mechanism in Cyg X- 
3, the space of possible ma-md combinations is still quite 
large, whereas RLOF as the mass-transfer mechanism leaves 
only few possible combinations of ma and toj- However, re- 
sults of population synthesis (see below) show that under the 
assumptions leading to the formation of He-starH-BH (HeH-BH) 
systems, only a minor part of the "allowed" area is populated. 
Fig.|2]shows that the observed mass-loss rate from the system 
and the observed P/P are hard to reconcile with the sugges- 
tions of donor masses of several IOA/q. 



sumed that the latter can accrete matter at M < MEdd and 
that the excess of the transferred mass is lost from the sys- 
tem, taking away the specific angular momentum of the ac- 
cretor. Prior to RLOF, mass loss by stellar wind was computed 
according to formulae Q. Accretion prior to RLOF was ne- 
glected. In the RLOF stage wind mass loss directly from the 
He star was neglected (see Fig.Q. The adopted ranges of ini- 
tial masses were 1.0Mq-4.1Mq for the He stars and 1.4Mq- 
S.OA/q for the compact objects. Computations were carried 
out using P. Eggleton's evolutionary co de (pr iv. comm. 2003, 
see also Eggleton & Kiseleva-Eggleton l2002l and references 
therein). A selection of the results are presented graphically in 
Fig. |3l The systems had the following combinations of com- 
ponent masses at the onset of mass transfer: — 3.0Mq, 
TOa = 5.OAf0 (Fig.lHa, b), Wd = 1.46Mq, = 1.4M0 
(Fig. El c, d), and ttia = I.OMq, = IAMq (Fig. El e, 
f). In all the computed systems the He stars started RLOF at 
P w 0.2 day. 

We find that the systems with He-star donors > S.OMq tra- 
verse the range of P/P observed for Cyg X-3 in ^ 10^ yr. We 
do note, however, that the typical value of P/P in the RLOF 
phase decreases with the mass of the donor (for a given mass of 
the compact star) and that the time spent close to the observed 
P/P range increases for lower-mass systems. A system that at 
the onset of mass transfer consists of a 1.46Af0 He star and a 
1.4Af0 CO. (Fig.|3lc, d) spends about 5 x IQ-^ yrs in the P/P 
range observed in Cyg X-3, and stays some 6 x 10"* yrs at P/P 
values less then twice those in the observed range. A system of 
a l.OA/0 He star and a 1.4Af0 c.o. stays within the observed 
range all throughout RLOF (Fig.|3]e, f). 



5. Population synthesis 

We carried out a population synthesis to determine the current 
number of He-nc.o. binaries in the Galaxy. The details of the 
population synthesis are bri efly d e scribe d in Appendix B. We 
us ed the a pproximations of Pols ( 1993) to the computations 
of IPaczvriski (,19_7_1) and .Habet& to estimate the core- 

He-burning times. In Figs. HI and [S] we plot the masses of the 
components and the orbital periods of the He-nc.o. systems that 
have He components in the core-He-burning stage. We find that 
there are currently 200 HeH-BH and ^ 540 He-nNS binaries 
in the Galaxy.^ 



4. Evolutionary caicuiations 

Analytical estimates for Roche-lobe-overflowing systems pre- 
sented above suggest that the "observed" M range of Cyg X- 
3 may be typical for moderate-mass (up to several Mq) He 
stars overflowing their Roche lobes after completion of core He 
burning. The results of our population synthesis also suggest 
that most He companions to compact objects are of moderate 
mass, see Sect.|5l 

To verify the inferences in the previous Section we car- 
ried out several evolutionary calculations of semidetached sys- 
tems consisting of He stars accompanied by a c.o. We as- 



5. 1. Wind- fed X-ray systems 

As noted in Sect. l3.2.'T1 we expect an accretion disk in the Cyg- 
X-3 system. The systems shown in Figs. |3 and [S] may form a 
disk through wind accretion if the wind matter carries enough 
angular moment um as realised first bv llUarionov & SunvaevI 
fcl975.) : see, e.g., iLiviol (17994) for later work on this subject. 



" Note that these numbers represent one possible random realisation 
of the model for a population of He+c.o. systems, so all numbers given 
are subject to Poisson noise. 
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Fig. 3. AI and P/P as function of Mhc for systems with Roche-lobe overflowing He-star donors and compact accretors. Masses 
(mdjTOa) at the onset of mass transfer are: (3.0, 5.0), (1.46, 1.4), (1.0, 1.4) Mq (from top to bottom). Roche-lobe overflow 
starts at P « 0.2 day. In all panels thick solid lines show results of computations. In panels for P/P thin solid lines show the 
limits of observed P/P in Cyg X-3. Dotted curves show rough estimates for M and P/P based on Eq. derived with the 
approximations to R and L at the terminal-age He main-sequence from Hurlev et al. (2000) . 
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5.1 .1 . Helium-star/black-hole binaries 

Following the derivation in ' Ergma & YungelsorJ ( ll99 Sl for sys- 
tems with Kerr BHs, we apply the disk-formation criterion 



P < O.2(MBH/M0)«r4 day, 



(4) 



where Mbh is the BH mass, and f looo is the magnitude of the 
radial velocity of the wind in the vicinity of the c.o. in units of 
1000 km s^^. It appears then that for wiooo — 1 there are cur- 
rently ~ 30 wind-fed disk-forming systems in the Galaxy, 5 of 
them with orbital periods similar to that of Cyg X-3. However, 
only 9 out of the 30 systems have A/hc ^ TMq donors and 
would be identified as WR stars. Only one of these systems has 
an orbital period close to that of Cyg X-3. The remaining four 
systems with P in the observed range have Mhc ~ (2 — A)Mq; 
according to Eq. (|2ji their wind mass-loss rate will be below 
~ lO~^Af0 yr~^ and they probably will not produce the WR 
phenomenon ^. There are 8 systems with Afne ^ ^M© and 
P > 10 hr that fulfil the disk-formation criterion. Note, how- 
ever, that criterion @ has a very steep dependence on fiooo; 
if we assume, e.g., that t'looo = 1-5, only 3 systems out of 
the 30 able to form disks remain, all of them with P < 0.15 
day and with A/ho ~ 3Af(n. Since k nown WR stars have 
wiooo ~ (0.7 - 5.0) jNu gis & LamerslEoO Oft. we claim that 
Cyg X-3 may well be the only wind-fed WRh-BH system in the 

Galaxy. 

As already n oted by llben etalJ ill 99,4 and 
lErgma & YungelsorJ \l99^, with vwoo ~ 1.5, Cyg X-3 
may have a wind-fed disk if A/bh ~ 5A/0. The latter value 
fits well into the model range of expected black-hole masses in 
wind-fed systems with disks and orbital periods close to that 
of Cyg X-3 (the star symbols in Fig.|4j. An additional reason 
why only a WR system with P < 0.25 days shows up as a 
WR X-ray binary may have to do with the velocity profile of 
the wind. In such a close system the wind will, at the orbit of 
the compact object, not yet have reached its terminal velocity 
(of for example 1500 km s^^), whereas in the wider systems it 
may have, such that no disk forms in the wider systems. 



5.1 .2. Helium-star/neutron-star binaries 

Figurelsjshows the population of He-nNS binaries, which can be 
divided into two subpopulations. The first and larger subpopu- 
lation consists of systems with P < 1 day and A/ho ^ 2A/0. 
The short period is due to their previous common-envelope 
(CE) phase; the large number of low-mass He-star systems is 
due to the initial-mass function and to the fact that low-mass 
He stars live much longer than higher-mass He stars. The other 
systems form through a "double spiral-in" in which two giants 
go through a CE phase, producing two He stars. This formation 
channel only occurs for nearly-equal-mass binaries, and since 



A possible reason for the absence of the WR phenomenon in 
low-mass He stars is that, for luminosity-to-mass ratios characteristic 
for them, one does not expect radial pulsations that may drive shock 
waves with a consequential increase of the g as density, n ecessary to 
produce radiative stellar winds (see, e.g.. iFadevev & Novikova.2003.) . 
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Fig. 4. Current population of core-He-burning HeH-BH bina- 
ries in the Galaxy. Upper panel — distribution in A/hc — P 
plane; middle panel — distribution in A/bh ~ P plane; lower 
panel — distribution in A/hc — A/bh plane. The dash-dotted 
vertical lines in the upper and lower panel show the lower- 
mass boundary for He stars identified as WR stars. Systems 
that satisfy the disk-formation criterion for wind-fed objects 
(Eq. 13 are marked by open circles; the subset of them with 
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Fig. 5. Current population of core-He-burning He-nNS binaries 
in the Galaxy. The dash-dotted vertical line shows the lower- 
mass boundary for He stars identified with WR stars. The four 
dashed lines show the critical periods below which according to 
criterion (|5} disk formation through wind accretion is possible, 
if "yiooo = 1^ 1-5, 2, and 3 (highest to lowest), respectively. 



the most massive He star collapses into a c.o. there is a lower 
limit to the mass of the other He star (see iBrownl 1 995l) . 

As shown by llUarionov & SunvaevI ll 19751) and 
lErgma & YungelsoiJ ( Il998l) an accretion disk will form 
in a HeH-NS binary if the rotation period of the NS is longer 
than the equilibrium period that was established during the 
CE episode that accompanied the formation of the He star 
Assuming that the wind-mass-loss rate may be described by a 
formula M = fcMg^ (e.g.|Nelemans & van den Heuvel 2001 ) 
one derives as the disk-formation criterion 

P < 1.5 • X^k^ '^^M^lM^^-^M^l^'^v^^QQ hr. (5) 

[This is equivalent with Eq. (9) oflEr gma & Yungelso n'('l99{4. 
but with a more general expression for M\. We assume fc = 
1.38 X 10"^, a = 2.87 after Nelemans & van den Heuvel 
(I2OOII) . which gives practically the same M values as Eq. (|2j- 
Limiting periods for uiooo = 1. 1-5, 2, and 3 are plotted in 
Fig. |5] (dashes). We find ~ 5 systems that may be identified 
as WRh-NS systems with accretion disks. However, the very 
steep dependence on the wind velocity may reduce this num- 
ber to 0. Thus, as already noted by Ersma & Yungelson ( 1998), 
the low angular momentum of WR-star winds may completely 
preclude the formation of wind-fed Cyg-X-3-like systems with 
NS companions. 

One should note, however, that the mass-loss rates for He 
stars quoted above were derived from observational data on 
WR stars. Data on mass-loss rates for lower-mass He stars are 
not available. It is therefore not clear whether M may be ex- 
trapolated to below ^ IMq. Hence, the validity of criterion 
(|5j below this mass is uncertain. 



5.2. X-ray systems powered by RLOF 

Another possibility for Cyg X-3 is that the system contains a 
He star of < 7Mq which transfers mass in the case BB of 
RLOF. Then the WR spectrum arises in re-ejected matter Also 
in the RLOF case we can estimate the number of systems we 
should currently observe as WR X-ray binaries. As mentioned 
in Sect. 0]the RLOF systems that can provide a P/P close 
to the observed range must initially have had a low-mass (< 
1.5 Mq) He-star donor. From the population-synthesis calcu- 
lations we expect ~ 500 such systems in the Galaxy at any 
time. Typically these systems live ~ 1.5 x 10^ yrs, of which 
they spend ^ 10^ yrs in the phase of Roche-lobe overflow. We 
thus expect only of order 1 such system in the Galaxy in the 
phase of Roche-lobe overflow at any time. 

As indicated by Fig.|5] the bulk of these systems has P < 
10 hr. The typical mass-transfer rates in these systems are in 
the range (1 — 3) x lO~^Af0 yr^^, and most of the transferred 
mass will be lost from the system through re-ejection. These 
rates are consistent with the lowest observational estimates of 
a'/ for Cyg X-3. 

6. Discussion 

6.1. The "missing" He+c.o. binaries 

Our population synthesis shows that there are several core- 
He-burning HeH-NS binaries and a few dozen core-He-burning 
HeH-BH binaries with He-star masses > 7Mq in our Galaxy. 
If we assume that all matter that passes through the so-called 
accretion radius 7-a = 2Gma./c^ (G the gravitational con- 
stant, c the speed of light) is accreted by the c.o., and that 
the gravitational potential energy of the accreted matter is con- 
verted into luminosity, all wind-accreting He-star binaries with 
^^Hc > 3Mq in our model population will have an intrinsic 
luminosity > 10"^® erg s^^ and should be observable as He-star 
X-ray sources. In this we also assume that Eq. (|2} holds down 
to low-mass He stars. 

Apart from Cyg X-3 a few WRh-c.o. candidates are re- 
ported, e.g., HD 197 406AVR 148 (Marchen ko et al. 1993), 
HP 191765AVR 134 jMOTel et alJll999l) . HP 104994AVR 46 
('Marchen ko et alJbOOoT However, it is still unclear whether 
the companions to these WR stars really have a relativistic na- 
ture; the systems lack the X-ray luminosity expected in such 
case. A low Lx can be reconsiled with, e.g., a spinning pulsar, 
that deflects the flow. 

The fact that we do not observe several tens of He-star X- 
ray binaries in the Galaxy may be due to self-absorption of the 
X-ray photons by the wind of the donor It turns out that for 
the binaries in our population-synthesis sample the minimum 
column density between the c.o. and Earth due to the He-star 
wind depends mainly on A/ho- This column density is > 10 
g cm^'^ for all sources with A/ho > 4:Mq, rendering these 
sources unobservable in X-rays at energies below 20 keV. We 
do note that INTEGRAL has discovered several sour ces at en- 
ergies > 20 keV, of which 40 are still unidentified dUbertinil 
I2OO5I) . These hard sources might well be the missing He-nc.o. 
binaries. 
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The derivation for the column density also applies to 
RLOF-accreting sources that spherically symmetrically throw 
out overflowing matter in excess to the Eddington rate. We 
saw in Sect. |3] that the mass-transfer rate for RLOF systems 
is > 5 X IO^^Mq yr~^, well above the Eddington rate for a 
solar-mass c.o. (MEdd ~ 8 x lO^^Af© yr'^ for a 1.4 Af© NS 
accreting pure He) and good enough for a minimum column 
density of ^ 10^ g cm~^. This may support the suggestion of 
a model for Cyg X-3 in which the excess matter is thrown out 
of the system equatorially instead of spherically symmetrically 
(Sect.|3} together with a very low inclination for the system. 

6.2. IC10X-1 andSS 433, two other Wolf-Rayet 
X-ray binaries? 

6.2.1. IC10X-1 

iBauer & Bran dt* (2004") and lClark & CrowtheJ ( |2004 find that 
the luminous X-ray source IC 10 X-1 [io.1_2.5keV = (2 — 4) x 
10^8 erg s^^] in the starburst galaxy ICIO is spatially coin- 
cident with W NE star [MAC92] 17-A (notation adopted from 
ICrowther etal...2003;) . Assuming [MAC92] 17-A to be the most 
probab le optical counterpart of ICIO X-1, Iciark & CrowtheJ 
1I2OO4 fit a model with a stellar temperature of 85 000 K, 
logCL/Lg) = 6.05, Af = 4 X 10'^ Mq yr~\ and a terminal 
wind velocity of 1750 km s^^ to the observed He II A4686 
and N V AA4603-20 emission. They infer a mass for the 
WR star of ^ 35 Af©, using the_WR mass - luminosity re- 
lations o f ISchaerer & Maedeg (Il992h . Allowing for clumping, 
IClark & Crowthen find that M is equivalent to a homogeneous 
ma ss-loss rate of ^ 1 x IG ^^Af© yr^^. 

iBauer & Brandt! ( |2004 note fliat ICIO X-1 is quite simi- 
lar to Cyg X-3 in terms of X-ray luminosity, spectrum, and 
variability. Thus, ICIO X-1 may be the first extragalactic ex- 
ample of a short-living WR X-ray binary similar to Cyg X-3. 
The identification of the optical counterpart to ICIO XI with a 
massive WR star suggests that the system is wind-fed.'* 

Note, however, that in the field of ICIO X-1 there are three 
other candidate optical coun terparts, with O- o r B-sp ectral 
types, to the X-ray source. IClark & CrowtheJ (l2004l) sug- 
gest that their wind mass loss is insufficient to explain the 
observed X-ray lu minosity with wind accretion. From this 
Iciark & Crowthe j argue that, if one of those three candidates 
is the optical counterpart to ICIO X-1, the system is Roche- 
lobe-overflow fed similar to LMC X-4 or LMC X-3. 

6.2.2. SS 433 

Also the variable X-ray and radio source SS 433 was sug- 
gested to be a WR X-ra y binary ( van den Heuvel et a l. 1980: 
iFuchs et alJ2002alE004l) . Van den Heuvel et al. suggested that 
SS 433 contains an evolved early-type star or a WR star, based 
on the nature of its stationary spectrum, the size of the emitting 

Note that the estimate of M quoted above is more than an order of 
magnitude lower than would be expected for such a massive WR star, 
even when keeping in mind the low metallicity of ICIO: Z/Zq ~ 0.25 
(e.g.. .Leaueux et al..l979i) . 



region, the necessary presence of a strong wind for the produc- 
tion of the IR emission, and the large outflow velocity of the 
wind. 

In one of the latest attempts to identi fy the optical coun- 
terpart of SS 433, Fuchs et al] ( l2002albl) compared its mid- 
IR spectrum to WR stars of the WN subtype. They found 
the spectrum of SS 433 to resemble that of WR 147, a 
W N8-HB0.5V binary with c olliding wind. Using the formulae 
of Wright & Barlow ( 1975) and taking wind clumping into ac- 
count, Fuchs et al. (2004) obtain M = (5-7) x IO^^A/q yr~\ 
compatible to a strong late-WN wind. 

IFuchs et alJ ( 1200 4) propose that the material suiTounding 
the c.o. forms a thick torus or envelope around it rather than 
a classical thin accretion disk. They argue that the material is 
ionised by X-rays emitted from the vicinity of the c.o. and ex- 
pelled by radiation pressure which results in the imitation of a 
WR star. As mentioned above, this model needs further elab- 
oration, especially the formation of the WR spectrum and the 
self-absorption of the X-rays. 

On the other hand, lKingetalJ(l200Cl) suggested that SS 433 
is a mass-transferring system in which the formation of a CE 
may be avoided if radiation pressure expels the transfeiTed mat- 
ter in excess of the Eddington rate, i.e., a re-ejection model with 
a hydrogen-rich donor. For this model the donor mass must be 
in the range of (4 — 12)Af0. The model received support by 
the discovery of A-type super-giant features i n the spectrum of 
SS 43 3 observable at certain orbital phases (iGies et aljEooi 
iHillwig et aL.2004.) . Estimated masses of the comp onents are 
10.9 ± 3J_and_2^ OTMo, fitting the lKing et all model. As 
noted by Fuchs et al." (2004), if the results of Gies et al.' ('2QQ^ 
and H illwig et ak ( .2004) are confirmed, one needs to resolve an 
apparent contradiction of simultaneous presence of A-star and 
WR-star features in the spectrum. 

If the presence of an A-type star is not confirmed, however, 
it may appear that SS 433 really is a WR X-ray binary, the 
second known in the Galaxy after Cyg X-3. 

7. Conclusions 

We find in Sect. [S] that in principle there are two possible 
He-star binary configurations which may explain the observed 
P/P and M of Cyg X-3, and for which population-synthesis 
calculations combined with disk-formation criteria predict the 
existence of ~ 1 such system in the Galaxy. 

The first possibility is a system consisting of a massive (> 
7Mq) helium (i.e., WR) star and a BH around which a disk is 
formed through wind accretion. Population synthesis predicts 
that at any time ~ 1 such system with an orbital period similar 
to that of Cyg X-3 is present in the Galaxy, provided that the 
wind velocity near the orbit of the c.o. is < 1000 km s~^. In 
this case the system will have a lifetime of several times 10^ 
yrs (the lifetime of the He star) and the secular orbital-period 
increase is simply due to stellar-wind mass loss. 

The second possibility is a system consisting of a He star 
with a mass < I.5A/0 and a NS, which is powered by mass 
transfer due to RLOF, at a rate in the range (1 — 3) x 10^^ Af© 
yr~^. Population synthesis predicts that also ~ 1 such system 
with F < 10 hrs may be present in the Galaxy at any time. In 
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this case the system will have a lifetime of order 10^ yrs and the 
secular orbital-period increase is due to the combined effects of 
the mass transfer and subsequent mass loss — at a rate close to 
the transfer rate — from the accretion disk around the NS. 

In view of the population-synthesis results, we deem 
both configurations equally likely for Cyg X-3. We note that, 
though the first of the solutions implies the presence of a "real" 
WR star in the system, its mass is probably not extremely 
high, < 12 Mq. Thus, both solutions are consistent with the 
conclusion of a relatively moderate mass for the companion in 
Cyg X-3, which follows from the identification of the emission 
regiori with the vicinity of the compact object jpender et alJ 
ll999HStark & Saial2003h . 

We now speculate on the fate of these configurations. 

In the "wind" case, if the WR star loses sufficient mass, 
it might terminate as a NS, such that a system consisting of 
a BH plus a young radio pulsar would emerge. In view of 
the likely mass of ~ 5Mq for the BH (a value well within 

the range of BH masses in co nfirmed BH binaries, see, e.g., 

I " i |i 
McClintock & Remillar d 2004), disruption of the system in the 

supernova explosion seems unlikely. Alternatively, the WR star 
might collapse to a black hole, producing a double-BH binary. 

The fate of the "RLOF" configuration is most likely a bi- 
nary consisting of a massive white dwarf (composed of CO or 
ONeMg) together with a recycled neutron star. Several such 
systems are known in our galaxy as fast-spinning binary radio 
pulsars w ith massive w hite-dwarf companions in circular orbits 
(see, e.g.. lStairsl2004l) . 
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APPENDIX A 

We derive below the dependence of P/P on the masses of 
the donor and the accretor md and ma, the mass-loss rate from 
the donor mj, the fraction of riid that goes in a direct wind a, 
and the fraction of rhd that is re-ejected after transfe r to the 
accreto r /3. Early similar derivations are given in e.g. iHuand 
( 1 1963b andEitukov & Yungelson ( 1971). 

The total orbital angular momentum of the system is 
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"P 



(ma 
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From Eq. (|6} we have 
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dJ 
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(6) 



(7) 



We then use Kepler's third law and definitions 

dma = —(1 — a — /3)dmd 



d(ma + md) = (a - 
(by definition, a + (3 < 1) to derive 
,dJ (1 — a — /3)dmd 



dP 
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■ /3)dmd 



drriA 



(8) 
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(10) 

The logarithmic derivative of the orbital angular momentum 
due to mass loss from the system is given by 



(11) 



d J aWa + (3m\ dra^ 
J (ma + md) mamd ' 

Inserting this into ( I10> we obtain 

dP dmd 

P mi^mdimi, + nid) 

X [(3ma — 2mamd — 3md)Q; — 2mamdP + 3md — 3m^] 

(12) 

We note that our Eq.^]is not consistent with the expres- 
sions for variations of a and P in a similar model given by 
ISoberman et al ] (ll997l their Eqs. 30 and 31). The correct equa- 
tions in their notation are: 

=2«,-l) + (l-2^„) " 



dlnq 



d\nP 
d\nq 



+ (3 + 2^„)e 
- 1) + (1 - 3^„ 
+ (5 + 3'^Oe 
APPENDIX B 



q 



l + eq 
''l + q 

q 

1 + eq' 



(13) 



(14) 



We used the iNelemans et al ] ll2004 realisation of the 
SeBa program dPortegies Zwart & Verbunll Il996l) to caii-y 
out our population-synthesis calculations. The effects of stel- 
lar wind and c.o. formation are shown in Fig. |6l The 
most important assumptions about the evolution of mas- 
si ve stars in binaries relevant for this paper ar e as detailed 
Portegies Zwart & Yungelson (1998. 1999): Nelemans et alJ 



m 



il200l[) with the following two differences: 

(1) WR stellar winds are modelled by the fit to observed 
WR m ass-loss rates derived by .Nelemans & van den Heuveil 
ll200ll) . 

(2) It is assumed that stars with carbon-oxygen cores more mas- 
sive than 6.5 Mq collapse into a BH (for single stars this means 
stars with an initial mass above 25 Mq, in close binaries stars 
with an initial mass above 28 Mq, see Fig.|6j. The mass of 
the resulting BH is 6 5 per cent of the mas s of the exploding 
object (based on, e.g.. lNelemans et alJl999l) . The lowest-mass 
BH (for an exploding naked carbon-oxygen core) thus has a 
mass of 4.2 AIq. 
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Fig. 6. Summary of the evolution of the mass of single stars 
or components of wide binaries (top panel) and of stars in bi- 
naries experiencing RLOF (bottom panel). As a function of 
initial mass the lines give respectively: the mass at the end of 
the main-sequence (Aftams, solid line), the initial mass of the 
He-core (Mcorc,i, dotted line), the mass just before the super- 
nova explosion or the formation of the white dwarf (Afpre-SN, 
dashed line) and the mass of the final c.o. (A/comp.obj, black- 
white-black line). The bottom panel shows the masses for a 
primary that loses its hydrogen envelope soon after the end of 
the main sequence, before the He-core burning starts. 



The bottom panel of Fig. |6l summarises the evolution of 
binaries in which RLOF occurs. For which period this hap- 
pens is a function of the mass of the primary. As can be seen 
from the figure we skip RLOF for any case-B period for initial 
masses > 80Mq. Note that RLOF may be omitted for much 
lowe r primary masses due to t he occurrence of an LBV phase 
(e.g.. lVanbeveren et akF lQQS?). However, a detailed discussion 
of this matter is beyond the scope of this paper. 
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